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Using the local density approximation (LDA) and a realistic phonon spectrum we determine the 
momentum and frequency dependence of a 2 F(]s.,uj) in YBa2Cu307 for the bonding, antibonding, 
and chain band. The resulting self-energy E is rather small near the Fermi surface. For instance, 
for the antibonding band the maximum of ReYl as a function of frequency is about 7 meV at the 
nodal point in the normal state and the ratio of bare and renormalized Fermi velocities 1.18. These 
values are a factor 3-5 too small compared to the experiment showing that only a small part of E 
can be attributed to phonons. Furthermore, the frequency dependence of the renormalization factor 
/7(k, uj) is smooth and has no anomalies at the observed kink frequencies which means that phonons 
cannot produce well-pronounced kinks in stoichiometric YBa2Cu3C>7, at least, within the LDA. 

PACS numbers: 71.10.Fd,71.38.-k,74.72.-h 



Angle resolved photoemission (ARPES) experiments 
in high-T c cuprates show that electrons near the Fermi 
surface interact strongly with bosonic excitations caus- 
ing a kink in their dispersion at rather well-defined 
energies [H 0, The most popular candidates for 

the bosonic excitations are spin fluctuations Q and 
phonons 0- Plausible arguments have been put forward 
for either of the two choices but no consensus has been 
reached up to now which of these two possibilities actu- 
ally is realized in the cuprates. Recently, an explanation 
of these effects entirely in terms of phonons has been 
proposed identifying the 40 and 70 meV kink energies 
with the buckling and the breathing phonon modes [5|. It 
was also argued that the electron-phonon (ep) coupling 
in these systems is strongly anisotropic with respect to 
momenta and phonon branches and that the coupling to 
the breathing and the buckling phonon modes takes place 
at different points in the Brillouin zonejBj], which would 
provide a natural explanation of several observations. 

Assuming that the LDA for stoichiometric YBa2Cu3C>7 
represents a reasonable first approximation for the 
electron-phonon coupling in doped cuprates it is the aim 
of this Letter to find out whether an explanation of the 
kink phenomenon entirely in terms of phonons is possi- 
ble. To this end we generalize previous calculations [6( for 
the Eliashberg function where both momenta were aver- 
aged over the Fermi surface to the case where the average 
is done only with respect to one momentum. Multiply- 
ing with 2/uj and considering a general electronic band 
v the resulting function A„(k, a;) can be used to study 
the anisotropy of the electron-phonon interaction as a 
function of the momentum and frequency using realistic 
phonon branches determined from first-principles. Our 
approach omits possible effects due to strong electronic 
correlations beyond the LDA, a problem which presently 
is unsettled. This shortcoming seems not to be very se- 
rious for our work and conclusions because several cal- 



culations suggest [3, la] that strong electronic correlations 
tend to decrease rather to increase the ep-interaction in 
q-integrated quantities such as the self-energy. 

The momentum and frequency dependent coupling 
function A 1/ (k, uj) is defined byQ 

2 

A„(k,w) = - V | ffj (ki/,k + q//)| 2 • 

UJ — 

<5(w-w q3 )<5(e k+qM ). (1) 

gj(\s.u, k+q/.t) denotes the amplitude for a transition from 
the electronic state with momentum k and band index 
v to the state with momentum k + q and band index 
/x creating (annihilating) a phonon with branch label j 
and momentum q (— q). ujqj and €k u are the eigenener- 
gies of the phonons and electrons, respectively, the latter 
are measured relative to the chemical potential. Addi- 
tional calculations suggest that A can be considered to 
be diagonal in the band indices. A direct evaluation of 
Eq. (fTJ) is difficult because high numerical costs restricted 
the calculation of the matrix elements g to rather coarse 
q-meshes. Large fluctuations due to a small number of 
terms in Eq.([T]) can be avoided, however, by introduc- 
ing Fermi-surface harmonics (Toj . In this approach, the 
k-dependence of A !/ (k, uj) is expanded in a basis set of 
functions W formed by products of an even power of elec- 
tron velocities u a (k^) = de^/dk-a where a is a cartesian 
index. We explicitly write 

A„(k, uj) = £V(n,ki/)A„(n,w) ) (2) 

n 

W(n, fa/) = v 2 x n * (k^ 2n » (ki/)t£ n « (ki/), (3) 

with n = [n x , n y ,n z ). Defining the Fermi surface average 
of a variable X(k) by 

((X))^ = J2S(e^)X(k)/J2S(e^), (4) 

k k 
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we define 



A„(n,w) = ((W(n,ki/)A„(k ) w)» I/ . 



(5) 



Inserting Eq.([2]) into Eq.© one can solve for the coeffi- 
cients Ai/(n, w) and one obtains, 

A„(k,w) = £W(n,ki/)£S- 1 (n,m)A I ,(m,w) ) (6) 

n m 

with 

S„(n,m) = ((If(n,b)lf(m,k I /)))„. 

Calculation of the quantities A„(n,w) involves Fermi- 
surface averages of the squared matrix elements g 
weighted by the basis functions W, which can be eas- 
ily performed within LDA. Once A^n, a;) are known, 
A„ (k, oj) can be obtained from Eq.©. Typically 10 func- 
tions W(n, ki/) were used in the following calculations. 
The sum over electron momenta in Eq.Q was carried 
out on a dense 36x36x4-mesh. For the momenta q we 
used a 12xl2x4-mesh. <5-functions in Eqs.(UJ and (0J in- 
volving electronic energies were approximated by a Gaus- 
sian with a width of 0.2 eV. Additional calculations us- 
ing more basis functions and a smaller width for the 
Gaussian showed that our results are reasonably stable 
against larger basis sets or a refined representation of the 
6- functions. 

Fig[TJ shows X u (k.,uj) integrated over the frequency 
from to 50 meV (left column) and from 50 to 80 meV 
(right column). The diagrams in the first row refer to 
the bonding, those in the second row to the antibonding, 
and those in the third row to the chain band. The lat- 
ter two are partly hybridized in the (tt,0) direction. The 
diagrams were calculated for k z = 0.125, but the corre- 
sponding ones for k z — 0.375 look very similar, i.e., the 
diagrams are practically independent of the momentum 
perpendicular to the planes. The thick and thin black 
lines correspond to lines with = and = ±0.2eV, 
respectively. As the above numerical procedure incorpo- 
rates only quantities in the vicinity of the Fermi surface, 
the calculated values for A„(k, u;) are reliable in the re- 
gion between the thin black lines whereas larger errors 
may be expected for momenta far away from the Fermi 
surface. 

The diagrams in the first row of FigfT] show that the 
coupling functions associated with the bonding band 
have approximately tetragonal symmetry but that they 
vary along the Fermi line: Their minima and maxima 
occur at the nodal and antinodal points, respectively, 
leading to a variation of about a factor 1.5 and 1.3 for 
phonons below and above 50 meV, respectively. Whereas 
a momentum-dependent A would be in line with the 
scenario of Ref.[5| there are serious discrepancies: The 
coupling to high-frequency phonons, for instance to the 
breathing and half-breathing modes, should show a max- 
imum in nodal direction whereas we find that it has 
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FIG. 1: Coupling function A I/ (k, u) integrated over to between 
and 50 meV (left diagrams) and u) larger than 50 meV (right 
diagrams). Shown are zy-cuts through the irreducible Bril- 
louin zone (0 < k x / y < n) for k z = 0.125. The upper plots 
refer to the bonding, the middle ones to the antibonding, 
and the lower ones to the chain band. The latter two are hy- 
bridized in the (n,0) direction. The grey areas represent states 
with energies far away (> 0.8 eV) from the Fermi energy. 



there a minimum similar as in the case of low-frequency 
modes. In other words, the coupling to high- and low- 
frequency modes at the nodal and antinodal points is 
very similar which means that electronic self-energies due 
to phonons at these two points will not select predom- 
inantly low- or predominantly high-frequency phonons. 
Furthermore, the variation of A along the Fermi line is 
rather small in our case, i.e., A cannot be considered to 
be very anisotropic. The diagrams in the second row of 
FigfT] describe A associated with the antibonding band, 
except for a part close to (ir, 0), where an exchange of 
character with the chain band occurs (lower diagrams of 
FigfT]). However, the antibonding states in the antinodal 
direction along k y w 7r and in the nodal direction are not 
affected by this hybridization. The momentum depen- 
dence of A associated with the antibonding band resem- 
bles that of the bonding band. The coupling of electrons 
near the Fermi line to low- and to high-frequency phonons 
is again rather similar as a function of k, the absolute val- 
ues for A are roughly larger by a factor 3 for low- than 
for high-frequency phonons, similar as in the case of the 
bonding band. The same trend is found for the chain 
band, with an anisotropy of the order of two or smaller. 
Using a 4x4x2-mesh for q yields qualitatively similar re- 
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FIG. 2: (Color online) Momentum and frequency resolved 
coupling function A(k, uf) in 1/meV of the antibonding band 
near the Fermi surface in the nodal (upper diagram) and the 
antinodal (lower diagram) direction. Solid and dashed lines 
correspond to k z — 0.375 and 0.125, respectively. 



suits for A as the 12xl2x4-mesh. Quantitatively, the dif- 
ference between bond and antibonding bands as well as 
their anisotropics are more pronounced in this case, in 
particular, for the bonding band. The present inversion 
procedure leads in this case also to unphysical small neg- 
ative values for A in a small region away from the Fermi 
line for the bonding band which does not occur for the 
employed finer 12xl2x4-mesh. 

The following results refer to the antibonding band, 
thus we drop from now on the band index v. Figf5] 
shows frequency-resolved coupling functions A(k, u>) for 
momenta near the Fermi line in nodal (upper diagram) 
and antinodal (lower diagram) directions. The solid lines 
correspond to k z — 0.375, the dashed lines to k z — 0.125. 
The differences between solid and dotted lines are very 
small reflecting the very small /c z -dependence of A. The 
spectra are very broad and extend over the entire one- 
phonon range. It is evident that models which use only 
one or a few selected phonons are not able to represent 
adequately the phonon spectrum. 

The imaginary part of the retarded electronic self- 
energy due to phonons can be written as [9] 



7m£(k, uj) 



duu\(k, u) 
2b{u) + f(u- uj) + J{u + uj) 



Here, b and / denote the Bose and the Fermi function, 
respectively, and we assumed the self-energy to be di- 
agonal in the band index. The corresponding real part 
of £ was obtained by a Kramers-Kronig transformation 
using a cutoff of 1 cV. FigJ3] shows the calculated self- 
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FIG. 3: (Color online) Real and imaginary parts of the self- 
energy E of the antibonding band near the Fermi points along 
the nodal and antinodal directions at the temperature T=l 
meV. 



energy near the nodal and antinodal Fermi points using 
the coupling functions of Figj^l The negative imaginary 
part of £ increases monotonically over the whole phonon 
spectrum and then approaches a constant. The nega- 
tive real part increases first, passes then through a maxi- 
mum roughly in the middle of the phonon spectrum and 
then decreases monotonically. Both the real and imag- 
inary parts of £ exhibit small oscillations as a function 
of frequency which reflect the many maxima in A(k, uj) 
and the complexity of the phonon spectrum. At elevated 
temperatures (for instance, at T=9 meV) these oscilla- 
tions vanish, the maximum in ReT. becomes somewhat 
smaller and JmS finite at zero frequency but otherwise 
there is no substantial change in S. FigJ3]also illustrates 
the small anisotropy of S: In the antinodal direction S 
is only slightly larger compared to the nodal direction 
while the position of the maximum of ReT, is approxi- 
mately the same. For the nodal direction experimental 
values for the maximum of i?eS are 20-40 meV[l, |Tl} 
whereas we have only about 7 meV at T=9 meV. 

Fig 0] shows the frequency dependence of the renormal- 
ization function Z, defined by Z(k, ui) — 1— i?eS(k, uj)/lo, 
near the nodal and antinodal Fermi points. At zero fre- 
quency Z is equal to 1 + A(k) where A(k) is the frequency- 
integrated, dimensionless coupling constant depending 
still on the momentum. This coupling constant is about 
0.18 and 0.22 at the nodal and antinodal points, respec- 
tively, illustrating again the smallness of the electron- 
phonon coupling. With increasing frequency Z increases 
first slightly, passes through a maximum and then ap- 
proaches unity beyond the end of the one-phonon spec- 
trum. Superimposed on this general behavior are many 
small oscillations which reflect the numerous density 
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ing band; b) It is important to keep a realistic phonon 
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FIG. 4: Renormalization function Z(k,aj) of the antibonding 
band near the Fermi points in nodal and antinodal directions 
for k z = 0.125 at the temperature T=l meV. 



FIG. 5: (Color online) heft: Electronic dispersion; right: 
spectral function of the antibonding band for different mo- 
menta k x = k y along the nodal direction using T=l meV and 
k z = 0.125. 



peaks in Fig[2] and which vanish already at T=9 meV. 
If, as in our case, the imaginary part of £ is rather 
small, the dispersion of the electrons is given by the 
equation uj = ek/Z(k, uj) which for small E's has the 
solution uj = ek/^(k, £k). A kink in the dispersion is 
generated if Z changes from a large to a small value in 
a narrow frequency interval. This happens, for instance, 
in a model with one single optical phonon, see Fig. 2 of 
Ref. [12j. However, for our realistic phonon spectrum and 
coupling functions this does not occur according to Fig0J 
so there will be no well-pronounced kink feature in the 
electronic dispersion. There seems to be a general con- 
sensus that ARPES data for the nodal and antinodal 
directions yield kink energies of about 70 and 40 meV, re- 
spectively. It is clear that these energies do not represent 
any characteristic energies in the frequency dependence 
of Z in Figg] 

Finally, the left diagram in Figj5] illustrates the weak- 
ness and smoothness of the renormalization of the disper- 
sion of the antibonding band in the nodal direction and 
also its weak dependence on temperature. The ratio of 
asymptotic slopes at small and large frequencies is about 
1.18 which is much smaller than typical experimental val- 
ues of about 2 in optimally doped cuprates [2I I3L [l3| . The 
right diagram in Fig[5] shows spectral functions of the 
antibonding band along the nodal direction for five dif- 
ferent momenta k x — k y . Unlike in the case of one op- 
tical phonon coupled strongly to electrons, discussed in 
Ref. [l2j , the solid lines show always only one peak which 
departs and then returns very smoothly from the nonin- 
teracting peak with decreasing energy. 

Based on the above LDA results our main conclusions 
are: a) The coupling function A(k, u>) shows as a func- 
tion of k a minimum at the nodal and a maximum at 
the antinodal Fermi point rather independently of the 
involved phonon modes. The resulting anisotropy in k is 
rather moderate both for the bonding and the antibond- 



spectrum and not only one or two phonon branches be- 
cause the w-dependence of A is mainly determined by 
phonon densities where, for instance, the breathing or 
the buckling mode do not play any distinguished role; c) 
The frequency- integrated coupling function A(k) varies 
between 0.18 and 0.22 along the Fermi surface. Due to 
these rather small values the phonon-induced electronic 
self-energies along the nodal direction are too small by 
about a factor 3-5 to account for the self-energies de- 
duced from ARPES experiments and for the observed 
slope changes in the electronic dispersion, at least, within 
the LDA. 
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nition of A in a multi-band superconductor. 
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